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ABSTRACT 

The origin and possible universality of the stellar initial mass function (IMF) is a major issue in astrophysics. One of the main objectives of the 
Herschel Gould belt survey is to clarify the link between the prestellar core mass function (CMF) and the IMF. We present and discuss the core 
mass function derived from Herschel data for the large population of prestellar cores discovered with SPIRE and PACS in the Aquila Rift cloud 
complex at d ~ 260 pc. We detect a total of 541 starless cores in the entire ~ 1 1 deg 2 area of the field imaged at 70-500 /im with SPIRE/PACS. 
Most of these cores appear to be gravitationally bound, and thus prestellar in nature. Our Herschel results confirm that the shape of the prestellar 
CMF resembles the stellar IMF, with much higher quality statistics than earlier submillimeter continuum ground-based surveys. 

Key words. ISM: individual objects (Aquila Rift complex) - Stars: formation 



1. Introduction 

The question of the origin and possible universality of the IMF, 
which is crucial for both star formation and galactic evolution, 
remains a major open problem in astrophysics. The Herschel 
Space Observatory (Pilbratt et al. 2010), equipped with its 
two imaging cameras SPIRE (Griffin et al. 2010) and PACS 
(Poglitsch et al. 2010), provides a unique tool to address this fun- 
damental issue. This problem and other key questions about the 
earliest phases of star formation and evolution are central to the 
scientific motivation for the Herschel Gould belt survey, which 
will image the nearby (d < 0.5 kpc) molecular cloud complexes 
of the Gould belt using SPIRE at 250-500 //m and PACS at 70- 
160 jum (cf. Andre & Saraceno 2005 and Andre et al. 2010). 

Starting with the 1 .2 mm continuum study of the Ophiuchus 
main cloud by Motte et al. (1998), several ground-based (sub)- 
millimeter dust continuum surveys of nearby, compact cluster- 
forming clouds such as p Ophiuchi, Serpens, and Orion B have 
uncovered 'complete' (but small) samples of prestellar cores 
whose associated core mass functions (CMF) resemble the stel- 
lar IMF (e.g., Motte et al. 1998; Johnstone et al. 2000; Enoch 
et al. 2006; Nutter & Ward-Thompson 2007; - see also Alves et 
al. 2007 and Andre et al. 2007). Albeit limited by small-number 
statistics at both ends of the CMF, these findings favor theoret- 
ical scenarios according to which the bulk of the IMF of solar- 
type stars is largely determined by pre-collapse cloud fragmen- 
tation, prior to the protostellar accretion phase (e.g., Hennebelle 
& Chabrier 2008). The problem of the origin of the IMF may 
thus largely reduce to a good understanding of the processes 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 

** Figures 3, 4, 5, and 6 are only available in electronic format via 
http://edpsciences.org 



responsible for the formation and evolution of prestellar cores 
within molecular clouds. Apart from limited statistics, another 
limitation of ground-based submillimeter continuum determi- 
nations of the CMF is that they rely on uncertain assumptions 
about the dust properties (temperature and emissivity). Herschel 
is now making it possible to dramatically improve the quality of 
the statistics and to reduce the core mass uncertainties by per- 
forming direct measurements of the dust temperatures. 

In this Letter, we use the initial results of the Gould belt sur- 
vey, obtained toward the Aquila Rift complex as part of the sci- 
ence demonstration phase (SDP) of Herschel, to present the first 
CMF derived from Herschel data. We discuss the robustness of 
the resulting CMF based on several tests and simulations. 

The Aquila molecular cloud complex is located at the high 
Galactic-longitude end of the Aquila Rift, in the neighbourhood 
of the Serpens star-forming region. Here we adopt a distance 
of 260 pc for the Aquila star-forming complex. For a detailed 
introduction of the region and a discussion of its distance, we 
refer the reader to Bontemps et al. (2010). 



2. Observations and data reduction 

Our Herschel observations of the Aquila Rift complex were 
taken on 24 October 2009 in the SPIRE/PACS parallel-mode. 
A common ~11 deg 2 area was covered by both SPIRE and 
PACS scan maps, with a scanning speed of 60"s _1 . The field 
was observed twice by performing cross-linked scans in two 
nearly orthogonal scan directions. The combination of the nom- 
inal and orthogonal coverages reduced the effects of l/f noise 
(see Sibthorpe et al. 2008). 

The PACS 70 pm and 160 pm data were reduced with HIPE 
(Ardila et al. 2010) version 3.0. Standard steps of the default 
pipeline were applied starting from the raw data (level 0). We 
used file version 1 flat-fielding and responsivity in the calibration 
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Fig. 1. (a) Column density map derived from SPIRE/PACS observations of Aquila. The subregion referred to as the main subfield in 
the text is marked by the white rectangle. The cyan triangles mark the positions of the 541 starless cores identified in the entire field. 
The locations of the HII regions W40 and MWC297/Sh62 are shown. The PDR region around W40 is framed by white polygon, 
while the dashed square outlines the small region shown in more detail in online Fig.|5^. (b) Same as (a) for the main subfield, with 
a total of 452 starless cores, marked by cyan triangles. 



tree, instead of the built-in version 3 of those. Therefore, the flux 
density scale was corrected with the corresponding responsivity 
correction factors: 1.78 at 70 pm, and 1.43 at 160 pm. Multi- 
resolution median transform (MMT) deglitching and second- 
order deglitching were also applied. Baselines were then sub- 
tracted by high-pass filtering, with a median filter width corre- 
sponding to the full length (180'-190') of the scan legs taken 
from HSPOT (Frayer et al. 2007). The baseline fits were per- 
formed on the 'level 1' data using interpolation for the masked 
bright sources not to over-subtract true sky emission. In this 
way, we removed stripes and preserved spatial scales and dif- 
fuse emission up to the size of the maps. The final PACS maps 
were created using the photProject task, which performs simple 
projection of the data cube on the map grid. 

Our SPIRE observations at 250, 350, and 500 pm were re- 
duced using HIPE version 2.0 and the pipeline scripts deliv- 
ered with this version. These scripts were modified, e.g., to 
include observations taken during the turnaround of the tele- 
scope. A median baseline was applied to the maps and the 
'naive' map-making method was used. Online Fig.|3]shows the 
500/350/250 pm SPIRE images. The PACS 160/70 pm images 
of the same field are shown in Bontemps et al. (2010). 

For SPIRE, the absolute calibration uncertanty is estimated 
to be -15% (Griffin et al. 2010), while for PACS the absolute 
flux accuracy is within 10% in the blue filter, and better than 20% 
in the red filter (Poglitsch et al. 2010). The in-flight calibration 
of the SPIRE instrument is described by Swinyard et al. (2010). 

Besides cross-correlating the SPIRE and PACS maps to 
test their relative astrometry, we compared the astrometry of 
the Herschel images with publicly-available Spitzer 8 pm and 
24 pm data, as well as high-positional accuracy (<1") 3 mm 
IRAM Plateau de Bure observations of a small field at the cen- 
ter of the Aquila main filament (Maury et al. in prep.). In this 
way, we corrected the Herschel images for small astrometric off- 
sets (< 6") remaining between the SPIRE and PACS maps, and 
achieved a final astrometric accuracy better than ~ 2". 



3. Results and analysis 

Compact sources were extracted from the SPIRE/PACS images 
using getsources, a multi-scale, multiwavelength source-finding 
algorithm briefly described in Men'shchikov et al. (2010). 
Several sets of extractions were obtained, including one for the 
entire field and one for the main subfield of the Aquila com- 
plex (see Fig.[TJ. At this early stage of the scientific exploitation 
of the Herschel survey, we only considered robust sources with 
significant {S IN > 7.5) detections in at least two SPIRE bands, 
especially since the significance of the sources depends slightly 
on the adopted set of extractions. 

For the Aquila main subfield (see Fig.[TJ, Spitzer 24 pm ob- 
servations were used in combination with PACS 70 pm data to 
distinguish between starless cores and young (proto)stellar ob- 
jects (YSOs). In this subfield, objects detected in emission above 
the 5cr level at 70 pm and/or 24 pm were classified as YSOs, 
while cores undetected in emission (or detected in absorption) at 
both 70 pm and 24 pm were classified as starless. This classifi- 
cation yielded 452 starless cores in the Aquila main subfield. 

Outside the main subfield, we had to rely only on our PACS 
70 pm data to distinguish between starless and protostellar cores. 
Based on the results obtained in the main subfield, we estimate 
that the lack of Spitzer 24 pm information leads only to a ~3% 
error in the classification. Altogether, we identified a total of 541 
starless cores and 201 embedded YSOs in the entire field. The 
YSOs include ~45-60 Class protostars depending on the se- 
lection criteria (see Bontemps et al. 2010). 

Based on our Herschel data, we constructed dust temperature 
(7d) and column density (Z) maps. To do this, we first smoothed 
all Herschel images to the 500 pm resolution (36.9") and repro- 
jected them to the same 6" pixel grid. Weighted spectral energy 
distributions (SEDs) were then constructed for all map pixels 
from the 5 observed SPIRE/PACS wavelengths. 

Assuming single-temperature dust emission, we fitted each 
SED by a grey-body function of the form I v - B v (T,i)(l - e~ Tr ), 
where 7 V is the observed surface brightness at frequency v, 
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Fig. 2. (a) Differential mass function (dN/dlogM) of the 452 starless cores in the Aquila main subfield, approximated with a lognor- 
mal fit (red curve). The error bars correspond to y/N statistical uncertrainties. The core sample is estimated to be complete down to 
~0.2-0.3 M Q , The lognormal fit has a peak at ~0.6 M and a standard deviation of ~0.42 in logioM. For comparison, the dash-dotted 
line shows the single-star IMF (e.g., Kroupa 2001), and the dashed curve corresponds to the unresolved system IMF by Chabrier 
(2005). The dotted line shows a power law of the form dAVdlogM °c M~ (1 6 , which is the typical mass distribution of low-density CO 
clumps (see Kramer et al. 1998) The high-mass end of the CMF is fitted by a power law (dN/dlogM oc M" 1 5±02 ), while the Salpeter 
IMF is dAVdlogM °c M~ L35 . (b) Same as (a) for the subset of 314 candidate prestellar cores identified in the Aquila main subfield. 
Here, the best-fit power law to the high-mass end of the bound cores CMF gives the same result (dN/dlogM <x M~ lA5±02 ), while the 
lognormal fit peaks at ~0.9 M and has a standard deviation of ~0.30. See text for discussion. 



t v - Kyi. is the dust optical depth, and k v is the dust opacity per 
unit (dust+gas) mass, which was approximated by the power law 
k v = 0.1 (v/1000 GHzf cm 2 /g (cf. Beckwith et al. 1990). The 
dust emissivity index/? was fixed to 2 (e.g., Hildebrand 1983). 

Each SED data point was weighted by 1/cr 2 , where the rms 
noise <x was estimated in an emission-free region of the map at 
the relevant wavelength, and the calibration uncertainties were 
also included. The two free parameters T^ and 2 were derived 
from the grey-body fit to the 5 Herschel data points for all pix- 
els for which the fit was successful. Map pixels for which the fit 
was unsuccessful or unreliable were assigned the median dust 
temperature of the successful fits. Likewise, the column den- 
sity along the line of sight to pixels with unreliable fits was 
estimated directly from the surface brightness measured at the 
longest wavelength with a reliable detection, assuming the me- 
dian dust temperature of the successful fits. 

A similar SED fitting procedure was employed to estimate 
the dust temperature, column density, and mass of each core. 
Here, the SEDs were constructed from the integrated flux densi- 
ties measured by getsources for the extracted sources. Ignoring 
the distance uncertainty (see discussion in Appendix A of Andre 
et al. 2010), the core mass uncertainty is typically a factor of ~ 2, 
mainly due to uncertainties in the dust opacity law (k v ). 

Monte Carlo simulations were carried out to estimate the 
completeness level of our SPIRE/PACS survey. We first con- 
structed clean maps of the background emission at all Herschel 
wavelengths by subtracting the emission of the compact sources 
identified with getsources from the SPIRE/PACS images of 
Aquila. We then inserted a population of ~ 700 model star- 
less cores and ~ 200 model protostars at random positions in 
the clean-background images to generate a full set of synthetic 
Herschel images of the Aquila region. The model cores were 
given a realistic mass distribution in the 0.01-10 M range and 
were assumed to follow a M oc R mass versus size relation. The 
emission from the synthetic cores was based on spherical dust ra- 
diative transfer models (Men'shchikov et al. in prep.). Compact 
source extraction of several sets of these synthetic skies was per- 



formed with getsources in the same way as for the observed im- 
ages. Based on these simulations, we estimate that our Herschel 
census of prestellar cores is 75% and 85% complete above a core 
mass of ~0.2 and ~0.3 M , respectively, in most of the field. 
Likewise, our census of embedded protostars is more than 90% 
complete down to L^ \ ~ 0.2 L . Our survey may however be 
less complete than these values in the high background region 
around the W40 PDR (see white polygon in Fig.QJ). 

4. Discussion and conclusions 

4.1. Prestellar nature of the Aquila starless cores 

In this paper, we follow the naming convention that a dense core 
is called prestellar if it is starless and gravitationally bound (cf. 
Andre et al. 2000, Di Francesco et al. 2007). In other words, 
prestellar cores represent the subset of starless cores that are 
most likely to form (proto)stars in the future. 

Strictly speaking, spectroscopic observations would be re- 
quired to derive virial masses for the cores and determine 
whether they are gravitationally bound or not. However, mil- 
limeter line observations in dense gas tracers such as N2H + 
show that thermal motions generally dominate over non- 
thermal motions in starless cores (e.g., Andre et al. 2007). 
Assuming that this is indeed the case for the Aquila cores ob- 
served here, we may use the critical Bonnor-Ebert (BE) mass, 
Mgg x 2.4 Rre a 2 /G, as a surrogate for the virial mass, where 
Rbe is the BE radius, a is the isothermal sound speed, and G 
is the gravitational constant. The critical BE mass may also be 
expressed as M™ ( w I.IS^jP^/ , where P e xt is the external 
pressure, which may be estimated as a function of the column 
density of the local background cloud, S c ;, as P ext w 0.88 G Z 2 , 
(McKee & Tan 2003). For each object, we derived two esti- 
mates of the BE mass: (1) MsE(Robs) as a function of the ob- 
served core radius assuming a gas temperature of 10 K, and (2) 
MbeC^ci) as a function of the local background column den- 
sity measured from both our source-subtracted Herschel im- 
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ages and the near-IR extinction map of Bontemps et al. (2010 
- see also Fig.|6j5). We then calculated BE mass ratios of q-be — 
max[MBE(Robs), MBE(^ci)]/M b s and selected candidate prestel- 
lar cores to be the subset of starless cores for which ctbe ^ 2. 
Based on this criterion, 314 (or ~ 69%) of the 452 starless cores 
of the main subfield and 341 (or ~ 63%) of the 541 starless cores 
of the entire field were found to be bound and classified as good 
candidate prestellar cores. These high fractions of bound objects 
are consistent with the locations of the Aquila starless cores in a 
mass versus size diagram (see online Fig.|4]for the main subfield 
and Fig. 4 of Andre et al. 2010 for the entire field). 

The self-gravitating character of most Herschel cores in 
Aquila is supported further by an independent analysis of their 
internal column density contrasts. We here define the internal 
column density contrast of a core as 'Z pea i c /l, core , where Y, pea k 
and 2 core are the peak and mean column densities of the core, 
respectively. Assuming optically thin dust emission at Herschel 
wavelengths and neglecting any dust temperature/opacity gradi- 
ent, the internal column density contrast can be estimated from 
the core intensity values in the same form, as /J 6 /7 V . The in- 
ternal column density contrast is closely related to the degree of 
central concentration of a core defined by Johnstone et al. (2000) 
as C — 1 - (S core fLpeak). Based on their radial intensity profiles 
(cf. online Fig.[5]S), the Aquila starless cores have a median in- 
ternal column density contrast ~ 4. For comparison, the internal 
column density contrast is larger than 3.6 for supercritical self- 
gravitating BE spheres, while it is only 1 .5 for a uniform-density, 
non-self-gravitating object (Johnstone et al. 2000). 

As a final check, we also estimated the typical column den- 
sity contrast of the Herschel cores over the local background. 
For a critically self-gravitating BE core, the mean column den- 
sity £gg l as 1.56 (Pext /G) 1 ^ 2 is expected to exceed the column 
density of the local background cloud E c ; by a factor of 1 .46, if 
Pext ~ 0.88 G I?, (McKee & Tan 2003). The candidate prestel- 
lar cores identified with Herschel have a median column density 
contrast of ~ 1.5 over the local background, which is consistent 
with the conclusion that they are self-gravitating. 



less cores, implying that our census of starless cores is less re- 
liable in this area. We thus constructed another version of the 
Aquila main subfield CMF, excluding the 83 cores identified to- 
ward the PDR region. The excluded region of unusually high 
infrared background emission (see in Fig.QJ) was defined using 
the dust temperature map shown in Bontemps et al. (2010). The 
high-mass end of this CMF can be fitted by a very similar power 
law to that of Fig. |2ji: dAVdlogM oc M~ ' ■ 5±a3 . 

Thanks to the large number of starless cores identified with 
Herschel in Aquila (541 cores in the entire field), we have been 
able to consider several core subsamples and construct a statis- 
tically meaningful CMF in each case. We confirm that the shape 
of the prestellar CMF resembles the stellar IMF, using data of far 
higher quality statistics than earlier submillimeter ground-based 
surveys and more accurate core masses. Based on simulations, 
we conclude that our mass distributions are robust and do not 
depend strongly on different sets of extracted sources. 

The column density maps shown in Fig.[T]and online Fig. 5 
illustrate the tight correlation between the spatial distribution 
of the prestellar cores and the filamentary structure revealed by 
the Herschel images (Men'shchikov et al. 2010). The difference 
between our two SDP fields (Aquila Rift and Polaris Flare), 
in terms of filamentary structure and core mass distribution, is 
highlighted and discussed in Andre et al. (2010). 
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4.2. Prestellar core mass function in Aquila 

Figure [2^ shows the mass distribution (CMF) of the 452 starless 
cores identified in the Aquila main subfield. A power-law fit to 
the high-mass end of this CMF gives dAVdlogM oc M~' - 5±0 - 2 for 
M core > 2 M , which is very close to the Salpeter power-law 
IMF (dAVdlogM oc M' 135 - Salpeter 1955). The CMF of the 
541 starless cores identified in the entire field, shown by Andre 
et al. (2010), has a very similar shape. 

The robustness of the derived CMF was tested by selecting 
several subsets of sources, based e.g., on their physical nature or 
locations in the SPIRE/PACS maps. In particular, we compared 
the mass spectrum of starless and prestellar cores. Figure^ was 
then compared with Fig. [2J3, which shows the differential mass 
function of 314 bound cores in the main subfield. We obtained 
the same best-fit power law fit to the high-mass end (dAVdlogM 
oc M^ 1 • 45±0 - 2 ) as for Fig. [2^. Only the low-mass end changed 
when selecting the subset of candidate prestellar cores. 

We also selected subsets of starless cores based on their lo- 
cations in the Aquila field. Along lines of sight to the HII region 
W40, the associated photon-dominated region (PDR, Shuping et 
al. 1999) is the source of very strong extended background emis- 
sion at all infrared wavelengths. The strong background emis- 
sion in the Spitzer 24 /mi and PACS 70 /mi images makes it 
more difficult to discriminate between YSOs and compact star- 
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Fig. 3. SPIRE 500 //m (top), 350 ^m (middle), and 250 //m (bottom) images of the Aquila SDP field. See details about data reduction 
and map making in Section 2. The corresponding PACS 160 //m and 70 fim images are shown in Bontemps et al. (2010). 
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Fig. 4. Mass versus size diagram comparing the locations of the 452 starless cores identified with Herschel in the Aquila main 
subfield to both models of critical isothermal Bonnor-Ebert spheres at T — 7 K and T — 20 K (black solid lines) and observed 
prestellar cores from the ground-based (sub)-mm continuum studies of p Ophiuchi and NGC2068/2071 by e.g., Motte et al. (1998). 
The masses of the Herschel cores were derived as explained in the text and their sizes measured at 250 /mi. The 314 candidate 
prestellar cores of the main subfield are shown as filled triangles, while the other starless cores are shown as open triangles. The 
mass-size correlation observed for diffuse CO clumps is also displayed (shaded yellow band - Elmegreen & Falgarone 1996). The 
typical (5cr) detection threshold of current ground-based (sub)mm (e.g., MAMBO, SCUBA) surveys at d — 150 pc is shown by the 
dashed curve. The estimated 5<x detection threshold of our SPIRE 250 /im observations is shown by the blue curve. The lower and 
upper red diagonal lines correspond to constant 10 21 crrT 2 and 10 22 crrT 2 mean column densities, respectively. The arrow indicates 
the global shift of the sources if a distance of 400 pc were adopted instead of 260 pc for the Aquila Rift complex (see also Appendix 
AofAndreetal. 2010). 
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Fig. 5. (a) Close-up column density image (selected from Fig.QJ) showing a close view of several starless cores identified with 
getsources. The black ellipses mark the major and minor FWHM sizes determined by getsources for these cores at X = 250pm. 
Two protostars are also shown by red stars, (b) Intensity profiles returned by getsources for the starless core marked by the arrow in 
Fig.|5^, at the five Herschel wavelengths. 
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Fig. 6. (a) Column density map of the Aquila entire field derived from Herschel data. The effective FWHM resolution is 36.9". 
Unlike in Fig.Q] a uniform offset N^ - 3.8 x 10 21 crrT 2 has been added in order to optimize the match with the near-IR extinction 
shown in (b). (Indeed, the Herschel mapping does not constrain the zero level of the background emission.) (b) Near-IR extinction 
map of the same field based on 2MASS data (see Bontemps et al. 2010) and expressed in units of column density, using the relation 
JV H , = 10 21 cirT 2 x A v . The resolution is 2'. 



